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THE MIXED PROBLEM IN LIPSCHITZ DOMAINS WITH 
GENERAL DECOMPOSITIONS OF THE BOUNDARY 

J.L. TAYLOR, K.A. OTT, AND R.M. BROWN 



Abstract. This paper continues the study of the mixed problem for the 
Laplacian. We consider a bounded Lipschitz domain Q C R", n > 2, with 
boundary that is decomposed as dQ = D U N , D and A'^ disjoint. We let A 
denote the boundary of D (relative to 9f!) and impose conditions on the di- 
mension and shape of A and the sets A^ and D. Under these geometric criteria, 
we show that there exists po > 1 depending on the domain Q such that for 
p in the interval (l,po)i the mixed problem with Neumann data in the space 
L^{N) and Dirichlet data in the Sobolev space W^'^{D) has a unique solution 
with the non-tangential maximal function of the gradient of the solution in 
L^idil). We also obtain results for p = 1 when the Dirichlet and Neumann 
data comes from Hardy spaces, and a result when the boundary data comes 
from weighted Sobolev spaces. 



1. Introduction 

In this paper we consider the mixed problem, or Zaremba's problem, for the 
Laplacian. Let be a bounded Lipschitz domain in R" and write d^l = D U N, 
where D is an open set of the boundary and = dfl \ D. We define the LP-mixed 
problem as the following boundary value problem 

Au = 0, infl 
u^fo, on D 



(1.1) 



= on N 

(Vu)* e LP{dn). 



Here, (Vw)* stands for the non-tangential maximal function of Vu. The normal 
derivative du/dv is defined as Vu • i/, where v is the outward unit normal vector 
defined a.e. on dft. Throughout the paper, all boundary values of u and du/dv are 
defined as non-tangential limits. See Section 2 for precise definitions. 

The study of the mixed problem is a natural continuation of the program of 
study of boundary value problems in Lipschitz domains which began over thirty 
years ago. Dahlberg \j\ treated the Dirichlet problem for the Laplacian in Lipschitz 
domains, while Jerison and Kenig jl5] treated the Neumann problem with boundary 
data in and the regularity problem with Dirichlet data having one derivative 
in L^. Verchota studied the regularity problem with Dirichlet data with one 
derivative in , and Dahlberg and Kenig studied the Neumann problem with 
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data [9]. The mixed boundary value problem in Lipschitz domains appears as an 
open problem in Kenig's CBMS lecture notes [TH Problem 3.2.15]. There is a 
large literature on boundary value problems in polyhedral domains and we do not 
attempt to summarize this work here. See the work of Baculiaet al. [1 for recent 
results for the mixed problem in polyhedral domains and additional references. 

Under mild restrictions on the boundary data we can use energy estimates to 
show that there exists a solution of the mixed problem with Vu in of the domain. 
Our goal in this paper is to obtain more regularity of the solution and, in particular, 
to show that Vu lies in LP{d^l). Brown j2j showed that the solution satisfies 
Vu G L?'{dQ) when the data /jv is in L'^{N) and fjj is in the Sobolev space 
W^'^{D) for a certain class of Lipschitz domains. Roughly speaking, his results 
hold when the Dirichlet and Neumann portions of the boundary meet at an angle 
strictly less than tt. In this same class of domains, Sykes and Brown ^31, obtain 
results for 1 < p < 2 and I. Mitrea and M. Mitrea [23] establish well-posedness in 
an essentially optimal range of function spaces. Lanzani, Capogna and Brown f20J 
establish results in two dimensional graph domains when the data comes from 
weighted L^-spaces and the Lipschitz constant is less than one. The aforementioned 
results rely on a variant of the Rellich identity. The Rellich identity cannot be used 
in the same way in general Lipschitz domains because it produces estimates in i^, 
and even in smooth domains simple examples show that we cannot expect to have 
solutions with gradient in L^{drt). 

Ott and Brown [26] establish conditions on 51, N, and D which ensure uniqueness 
of solutions of the i^-mixed problem and they also establish conditions on ft, N, 
D and and fo which ensure that solutions to the L^-mixed problem exist. All 
of this work is done under an additional geometric assumption on the boundary of 
D. More specifically, the authors address solvability of the mixed problem for the 
Laplacian in bounded Lipschitz domains under the assumption that the boundary 
between D and N (relative to dil) is locally given by a Lipschitz graph. Under 
these conditions on ft, N, and D, they prove that there exists po > 1 depending 
on the Lipschitz constant of the domain and on the dimension n, such that for p 
in the interval 1 < p < po, the L^'-mixed problem with Neumann data in L^^N) 
and Dirichlet data in the Sobolev space W^'P{D) has a solution and this solution 
is unique in the class of functions satisfying (Vm)* G Lp^OVl). In the case p — 
1, they prove results for the mixed problem with data from Hardy spaces. The 
novelty of this paper is to address existence and uniqueness of solutions of the L^- 
mixed problem under more general conditions on the decomposition of the boundary 
into sets N and D. Our proof relies on a technique of Shen [28] to use reverse 
Holder inequalities to establish existence of solutions to the L^'-mixed problem. 
This technique allows for an immediate extension to the mixed boundary value 
problem with data from weighted spaces. We carry out a study of the mixed 
problem in weighted spaces in Section 7. As one step of this argument, we need to 
consider the regularity problem with boundary data in weighted L^-spaces. 

The boundary between D and N is an important feature of the domain in the 
study of the mixed problem. Assume that 13 is a relatively open subset of and 
let A denote the boundary of D (relative to dVl). Before stating our assumptions on 
ri, TV, and D, we introduce the following notation. We will use 6{y) = dist(?/, A) to 
denote the distance from a point y to A. Let Br{x) = {y : \y ~ x\ < r} denote the 
standard bah in R" and let *r(a;) = Br{x)r\^l. For x G let /^r{x) = Br{x)C]d^l 
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denote a surface ball. We note that the term surface baU is not ideal since the 
"center" x may not he on the boundary. In addition, we wiU need to be careful in 
places because Ar{x) may not be a connected set. See Section 2 for other relevant 
definitions. 

Our assumptions on D, and D are stated here. We will obtain results only when 
the parameter e in ()1.3|) is small. See Section 2 for a definition of the constant tq. 

(1.2) H, C R",7i > 2, is a bounded Lipschitz domain of constant M. 

The set A is an Ahlfors (n — 2 + €)~regular set: There exists M > such that 

(1.3) •H"-^+'(A^(a;) n A) < Mr"-2+^ for aU a; e A, < r < ro, 

with e > 0. The notation H"^^^'^(i5) denotes the (n — 2 + e)-dimensional Hausdorff 
measure of a set E. Our third main assumption is that the set D satisfies the 
corkscrew condition relative to d^. There exists M > such that 

(1.4) for all x e A, < r < tq, there exists x G D such that jx — a;| < r 
and 6{x) > M~'^r . 

Several previously studied cases of the mixed problem fall under the conditions of 
assumptions (|1.2p , (|1.3[) , and (|1.4p . Venouziou and Verchota [32] establish a solution 
to the L^-mixed problem in polyhedral domains in R"^. In one particular case, 
they are able to solve the mixed boundary value problem in the pyramid in R'^, 
when Dirichlet and Neumann data are assigned to alternating faces. This example 
is not covered by the earlier work of Ott and Brown [55^ because at the apex of the 
pyramid, the boundary between D and N is not locally given by a Lipschitz graph. 
The pyramid example is covered by the results in this paper. Another example that 
is covered by this work, but not the earlier work of Ott and Brown f2^, is the case 
where the boundary of I? is a Koch snowflake of dimension slightly greater than 
n-2. 

We now state the main theorem of the paper. The definitions are given in Section 
m Since we do not assume that the Dirichlet set D C dfl is an extension domain 
for Sobolev spaces, note that we must assume that the Sobolev space W^''^{D) is 
defined by restricting elements in W^'P{d^) to D. See Section [2] for a discussion of 
the constants in the estimates of this theorem. 

Theorem 1.1. Letfl and D satisfy conditions ll.2\} . il.3\) . and There exists 

an exponent qo > 2, which depends on M and n, such that if A satisfies il.S]) . with 
< e < (go ^ 2)/(go — 1),. such that the following statements are true. 
If P ^ li L^-mixed problem has at most one solution. 

b) If fN lies in H^{N) and fu lies in H^'^{D), the L^-mixed problem has a 
solution u which satisfies the estimate 

ll(Vu)*||Li(ao) < C{\\fN\\m(N) + WfoWm^iD))- 

c) If Po = '?o((l — e)/(2 — e)) > 1, then for p in the interval (l,po) the following 
holds: If fN lies in LP{N) and fn lies in W^'P{D), then the -mixed problem has 
a solution u which satisfies 

||(Vu)*||ip(gf2) < C{\\fN\\LP(N) + ||./d||wi.p(D))- 

The proof of the main theorem will proceed as follows. We begin by recognizing 
that it suffices to prove Theorem 11.11 in the case where the Dirichlet data is zero. 
This is because non-tangential maximal function estimates for the gradient of the 
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solution to the Dirichlet problem are known. When the Dirichlet data comes from 
a Sobolev space, these estimates were treated for p = 2 by Jerison and Kenig [15] 
and for 1 < p < 2 by Verchota [331 IM] • The case of the Dirichlet problem with 
data from a Hardy space was treated by Dahlberg and Kenig [9, and by D. Mitrea 
in two dimensions |23l Theorem 3.6]. 

The first main result presented in the paper is existence of solutions of the mixed 
problem when the Neumann data is an atom for a Hardy space. The proof of this 
result is contained in Section 4. The key ingredient of the proof is an estimate of 
the Green function for the mixed problem, which is proved in Section 3. In Section 
5 we prove uniqueness of solutions to the L^-mixed problem, p > 1. Section 6 
contains the proof of the result and Section 7 contains the proof of the weighted 
result. 

Acknowledgment. The authors thank the referee for their helpful remarks. Part 
of this work was carried out while Russell Brown was visiting the Mathematical 
Sciences Research Institute in Berkeley, California, whose hospitality is gratefully 
acknowledged. 

2. Preliminaries 

Throughout the paper, we will work under the assumption (II. 2[) that is a 
bounded Lipschitz domain. A bounded, connected open set is called a Lipschitz 
domain with Lipschitz constant M if locally i7 is a domain which lies above the 
graph of a Lipschitz function. More precisely, for M > 0, x G and r > 0, define 
a coordinate cylinder Zr{x) to be Zr{x) — {y : \y' — x'\ < r, — Xn\ < (1 + M)r}. 
Use coordinates (x',x„) G R"^^ x R and assume that this coordinate system is a 
translation and rotation of the standard coordinates. Then f2 is a Lipschitz domain 
if for each x e dfl there exists a coordinate cylinder and a Lipschitz function 
(f) : R"^^ R with Lipschitz constant M such that 

QnZrix) = {{y',yn) ■.yn> Hy')}(^Zr{x), 

dnnZrix) = {{y',yn):yn = Hy')}(^Zri^)- 

Fix a covering of the boundary by coordinate cylinders {Zri{xi)}fLi such that each 
^loOrV i+M^ i^i) ^^^'^ ^ coordinate cylinder. Let rp = min{ri : i = 1, . . . ,L}. 

For a Lipschitz domain D, we define a decomposition of the boundary for the mixed 
problem, dil. = D U N as follows. Assume that D is a relatively open subset of dil, 
N = dft \ D, and let A be the boundary of D (relative to d^l). The assumptions 
on the decomposition of the boundary for the mixed problem are given in ()1.3|) 
and (|1.4|) . Recall that 5{y) — dist(j/. A) denotes the distance from a point y to the 
boundary of D. 

Many of our estimates will be of a local, scale invariant nature and will hold for 
r less than a multiple of rg , and with a constant that depends only on the constant 
M in assumptions (ll.2p - (|1.4p . e in (jl.3p . the dimension n, and any L^'-indices that 
appear in the estimate. We say that an estimate depends on the global character 
of the domain if it depends on the above and also on the collection of coordinate 
cylinders which cover dil and the constant in the coercivity condition (|3.2I) . The 
notation A ^ B will mean that c~^B < A < cB for some constant c depending 
only on M and n. 

We now prove several consequences of the conditions ()1.3|) and ()1.4|) that we will 
appeal to later in the paper. 
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Lemma 2.1. Let satisfy U.2]) and let r be such that < r < vq. If x ^ dfl 

satisfies S{x) > r\J\ + M^, then Ar(x) C N or Ar{x) C D. 

Proof. Let x G dfl and fix y G Ar{x). Since < r < rg, we may find a coordinate 
cylinder Z which contains Ar{x). Let (j) be the function whose graph gives dfl 
near Z. Since y £ Ar{x), we have that \y' — x'\ < r. Now define a function 
x' : [0,1] R."~\ x'{t) = {l-t)x' +ty'. Then -f{t) = {x' {t),(l){x' (t) )) gives a path 
contained in dfl that connects x and ?/ and has length less than ryl + M^. Since 
5{x) > r\/l + M^, and (5 is Lipschitz with constant one, we have that S{'-f{t)) > 
for < t < 1. Since 7(t) does not pass through A, we have that both x and y 
must lie in either D or N. The point y was arbitrarily chosen in Ar{x), therefore 

Ar{x) CD 01 Ar{x) C N. □ 

The following lemma is adapted from a result found in Lehrback [U pp. 254- 
255]. Here and throughout the paper, we use cr to denote surface measure. 

Lemma 2.2. Let il and D satisfy U.2\) and U.S]) and let r satisfy < r < rg. 
Then for each x Cz A and < t < 2r, 

(2.1) cr (A^(a;) n {y : 5{y) < t}) < Ct^-'r"~^+'. 

Proof. Fix X S A, < r < rg and t with < t < 2r. By a standard covering lemma, 
there exists a finite, disjoint collection of surface balls {At{yi)} with yi S Ar(x)n A, 
i = 1, 2, . . . , m, such that {Ar{x) n A) C Ui=i ^3t{yi)- Using that the collection of 
balls {At{yi)}^i is disjoint and the (n — 2 + e)-regularity of A given in ()1.3p . we 
have 



< C^H"-2+^(At(y,)nA) 



1=1 

< cn"-^+' {A3r{x) n A) 

This calculation yields the estimate m < mt = C(r/t)"^^+^ for t < 2r, where C 
depends on M and the dimension n. 
Now, 

a{Ar{x)n{y:d{y)<t}) < a (A4t(2/,)) 

i=l 

< Ct"-^mt 

which proves the Lemma. □ 

Lemma 2.3. Let and D satisfy \1.2\) and \1.4^ . There exists a constant c such 
that if x £ D and < r < tq, then there exists xd G D such that \x — xo \ < r and 
Acrixo) C D. Furthermore, 

(2.2) a{Arix) n D) > cr''-\ 

Proof. Let x E D and < r < rp. We will break up the proof into two cases. First 
suppose that r/2 < d{x), and let xd = x. Let Zrixn) be the coordinate cylinder 
centered at xd with radius r. Using Lemma |2. II we see that S{xd) > r/2 implies 
that Ac,r{x) C D for ci = l/{2Vl + M^). 



ri-2+e 

r 
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Now consider the case where 6{x) < r/2. Accordmg to ()1.4p and Lemma [2.11 
given i S A, there exists x £ D such that |x — a;| < r/2 and Ac^rix) C D 
with C2 = 1/(2Ma/1 + M2). Recall that S{x) < r/2, and choose i on A so that 
|a; — = d{x). By the remark above, there exists x with |5: — xj < r/2 and 
Ac-^r{x) C D. Let XD — X. Then 

|a; — xd| < |x — 2:| + ji — a;£)| 

< r. 

Thus, we obtain Acr{x) <Z D ii c < min(ci,C2). Conclusion (j2.2p is an immediate 
consequence, after perhaps making c smaller. □ 

The next two lemmas establish integrability of the distance function 5. 

Lemma 2.4. Let Q, and D satisfy U.^) . il.S]) . and If x £ dft and r satisfies 

< r < ro, then for —1 + e < s < oo, 

S{yy da « r"-i max(r, S{x)y. 

Proof. Fix X £ dfl, r in the interval (0, ro) and consider several cases: a) s > and 
S{x) < r/4, b) s > and S{x) > r/4, c) s < and S{x) < 4r, and d) s < and 
S{x) > 4r. 

Case a) Assume that s > 0, 6{x) < r/4. Since 6{x) < r/4 and S is Lipschitz with 
constant one, we have S{y) < 5r/4 for all y € Ar{x). As 9r2 is the boundary of a 
Lipschitz domain, we have a{Ar{x)) ~ r"^^. Thus when s > 0, the upper bound 
/a,(x) '^(2^)' ^ C'r"-i+'* follows easily 

To obtain a lower bound, we begin by finding x £ A such that |a: — ij = S{x) < 
r/4. Then the corkscrew condition for D (|1.4p gives a point x with |a; — i| < 
r/4 and S{x) > r/{AM). If y £ Ar/i^sM){x) we have that S{y) > r/(8M), and 
^r/{8Ai){x) C Ar{x) since |a; — < r/2. Thus the lower bound ^(y)" c^cr > 

[a ^ivY da > Cr"-i+^ follows. 

Case b) Assume that s > 0, S{x) > r/4. Observe that if y G Ar{x), then 
S{y) < S{x) + r < 5S{x) and the upper bound ^.^^ i^ly)" da < r"^^(5(a;)^ follows. 

To obtain the lower bound, we use that if 5{x) > r/4 and y £ Aj./s{x), then 
S{x) < 2S{y) and hence we have ^{y)'' da > ^^^^ <5(2/)* rfcr > cr"~^(5(a;)^. 

Case c) Assume that — 1 + e < s < 0, 6{x) < 4r. We divide the surface ball 
Ar {x) using level sets of the distance function and then use Lemma 12.21 to obtain 

oo „ 

siyyda = Yl / ^^yy^^'^ 

A,^{x) J A^{x)n{y:2-f'-'^r<S{y)<2-'^r} 

oo 

< C^r"^2+"(2"''r)"+i"' 

fc=0 

< Cr"~i+^ 

In the last inequality above we use the assumption that s > — 1 + e to sum the 
geometric series. 

To obtain the lower bound, we observe that if S{x) < 4r, then for y £ Ar{x), 
S{y) < 5{x)-^r < 5r. Since s < 0, it follows that 6{yy > {5ry and the lower bound 
Ia^x) ^(yy do- > Cr"-i+'' follows easily. 
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Case d) Assume that — 1 + e < s < 0, 5{x) > 4r. If y G Ar{x), then we have 
3<5(x)/4 < S{y) < 5S{x)/4. Thus, we obtam 6{yy da « r"-M(a;)^ 

The result of the Lemma follows easily from the four cases above. □ 

The integrability of the function S over interior balls is a straightforward adap- 
tation of the previous result. 

Lemma 2.5. Let fl and D satisfy hi. 2(1 . il.3\) . and 11-4% and let r satisfy < r < 
ro- Then for s ^ {—2 + e,co), 

[ 6{yYdyf^r"~^max{r,6{x)Y. 

Throughout this work, the main tool for estimating solutions will be the non- 
tangential maximal function. Fix a > and for x £ 957, the non-tangential ap- 
proach region is defined by 

T{x) ^ {y en:\x-y\<{l + a) dist(y, dVl)}. 

Given a function u defined on 57, the non-tangential maximal function is defined 
as 

u*{x) — sup X e dVl. 

yer{x) 

We will also utilize a truncated non-tangential approach region, 

Trix)=r{x)f]Br{x), 

and, respectively, a truncated non-tangential maximal function, 
u*{x)— sup X e 9f2. 

It is well known that for different values of a, the non-tangential maximal functions 
have comparable L^'-norms. Thus, we suppress the value of a in our notation. 

The restrictions of u and Vu to the boundary in (II. ip are understood as non- 
tangential limits. Precisely, for a function v defined on D, and x € dfl, 

v{x) — lim v{y), 

provided that the limit exists. It is well-known that for a Lipschitz domain and 
V a harmonic function in 57, the non-tangential limits exist at almost every point 
where the non-tangential maximal function is finite. 

We now recall the definitions of atoms and atomic Hardy spaces. A function a 
is an atom for 957 if suppa C Ar{x) for some x € 957, ||a||ioo(aQ) < l/a{Ar{x)) 
and ada = 0. In our treatment of the mixed problem, we will consider atoms 
for the subset N C 957. We say that d is an atom for if a is the restriction 
to A^ of an function a which is an atom for 957. The Hardy space H^{N), where 
A^ C 957, is the collection of functions / which can be represented as ^jO-j, 
where each aj is an atom for A^ and the coefficients Xj satisfy J^j \ < oo. In the 
case where A^ = 957 this definition gives the standard definition of the Hardy space 
-ff^(957). The Hardy-Sobolev space i?^'^(957) is defined as the set of functions 
with one derivative in if ^(957). More precisely, we say that a function A is an 
atom for H^'^{dil) if A is supported in a surface ball Ar{x) for some x S 957 and 
II VtA|jioo(aQ) < l/a{Ar{x)). If ?; is a smooth function defined in a neighborhood 
of 952 then the tangential gradient of v is defined as Vtw — Vu — (Vw • v)]^, where 
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V is the outward unit normal vector. Then A is an atom for H^'^{D) if A is the 
restriction to D of an atom A in H^ '^{dVl). The space H^'^{D) is the collection of 
all functions which can be represented as ^j^j where each Aj is an element of 
H^^^{D) and Y.j < oo. 

Finally, we define the Sobolev space W^'P{dft) to be the collection of functions 
in LP{dil) whose tangential gradient also lies in LP{dfl). 

3. Green function estimates and reverse Holder inequalities 

An important step in the proof of the main theorem is to show decay of the 
solution to the mixed problem with Neumann data an H^(N) atom as we move 
away from the support of the atom. This decay is encoded in estimates for the Green 
function for the mixed problem which are proved in this section. The argument 
that ensues only requires that be a Lipschitz domain and that D satisfies (|2.2p . 

When working near the boundary, we will want to assume that part of the bound- 
ary is flat. This can always be arranged in a Lipschitz domain by flattening the 
boundary with a change of coordinates. Since flattening the boundary will change 
the coefficients, we need to consider operators L with bounded and measurable coef- 
ficients. Assume that L = div^V, and assume that the coefficient matrix A is real, 
bounded, and measurable, satisfles A* = A, and satisfies the ellipticity condition 
that for every ^ e R^, there exists a A > such that 

A|^p<^^^<A-l|cl^ 

The optimal A for which the above condition holds is called the ellipticity constant 
for L. 

We now define a weak formulation of the mixed problem for solutions of diver- 
gence form operators whose gradients lie in L'^{fl). Our goal is to prove that under 
appropriate assumptions on the data, the weak solution will have a gradient in 
LP{dQ) for 1 < p < po, for some po > 1- For fc = 1,2, . . ., W'^--p{Q) denotes the 
Sobolev space of functions having k derivatives in LP{n). For D an open subset of 
the boundary, let W^'^{fl) be the closure in W^''^{i}) of functions in C^(R") whose 
support is disjoint from D. Let w]^^^'^{dfl) be the restrictions to dil of functions 
in W^^in) and define Wjj^^^ '^{dn) to be the dual of W^^'^'^id^). We assume that 

the Dirichlet data is zero and the Neumann data lies in the space Wj~i^^^'^ (dfl). 
Consider the mixed problem 

{div^Vu = 0, inn 
u = 0, on D 

AVu-v^Jn, on iV. 

We say that u is a weak solution of the problem (|3.ip if u G T4^^'^(f2) and 

To establish the existence of weak solutions to the mixed problem, we assume the 
following coercivity condition 

(3.2) \\u\\L-ii-i)<c\\Wu\\L2(n), u€W]j\n). 

Under this assumption, the existence and uniqueness of weak solutions to the 
boundary value problem (|3.ip are a consequence of the Lax-Milgram theorem. In 
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our applications, fl will be a connected, bounded Lipschitz domain and D will be 
an open subset of the boundary. These assumptions are sufficient to ensure that 
holds. 

We also need to define a weak solution of the mixed problem on a subset of ft. 
Let il' be an open subset of fl. We say that u is a weak solution to Lu — f in ft' 
with zero boundary data for the mixed problem on dfl' D dfl if m G ^dn'noi^') 
and for each test function (j> which lies in W^an'riCDun") (^')' ^® have 



Solutions of the mixed problem are bounded. If m is a solution of Lu = in the 
domain "^rix) — Br{x) n fi, x G il, and u has zero data for the mixed problem on 
dft n d'i>r{x), then there exists a constant C > such that 



This may be proved by the Moser iteration method [IS], for example. 

Finally, we give an estimate on the boundary Holder continuity of solutions of 
the mixed problem. For this estimate we consider domains 5'r(a;) = Br{x) with 
X ^ ft and r < dmt{x, dft), or 5'r(a;) = Br{x) n Vt with x G dft and r < vq. In the 
second case, we assume that dfl D Br{x) lies in a hyperplane. A study of Holder 
continuity of solutions of elliptic equations may be found in work of Stampacchia 
[29] . Stampacchia gives a general framework for studying Holder continuity of 
solutions using the method of De Giorgi TO^. This framework is applied to the 
mixed problem in the case where the boundary of D is the bi-Lipschitz image of a 
hyperplane of co-dimension 2. Our assumptions allow for more general subdivisions 
of the boundary. 

Theorem 3.1. Let x € dVl and assume that < r < tq. Let u be a weak solution of 
the mixed problem in ^'r(a;) with zero data for the mixed problem on d'^r{x) H d^. 
Then there exists an exponent /3 > such that 



\ J <i/^(x) 

The constant C and the exponent /3 depend only on the ellipticity constant A as well 
as M and n. 

The proof here will follow the method of de Giorgi [10] as given in the monograph 
of Ladyzhenskaya and Ural'tseva [19]. Fix \l//{(a;o) as above. We say that a bounded 
function u lies in the space B{'ifR{xQ),'-f) if for each ^s(2^) C 5'_r(xo), a £ (0, 1) and 
k as below we have 



Here, Ak^s = {y G '^six) : u{y) > k} with k an arbitrary real number if d'^six) H 
£> = and fc > if 5*,(a;) n £» 7^ 0. 

Lemma 3.2. Let x e dft and let r satisfy < r < ro- If u is a solution of the 
elliptic operator L in ^r(a^) with zero data for the mixed problem on d'^r{x) n dft, 
then u G 7) and 7 depends only on the ellipticity constant for L. 




(3.3) 
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Proof. Fix '^s{y) which is contained in 'i>r{x) and a E (0, 1). Let 77 be a smooth 
cutoff function that is supported in Bs{y), with = 1 on Bs-as{y) and satisfies 
|V?7| < Co/{as). Let u+ = niax(u,0) denote the positive part of a function u. If k 
is as in the definition of the space B, then (u — fc)+r/^ may be used as a test function 
in the weak formulation of divAVw — and thus we have 

/ Ty^AVu • V(-u - fc)+ dy = -2 / r]{u - fc)+ AVu • Vt^ dy. 

Using the symmetry and non-negativity of A and Young's inequahty, we obtain 

2/ ifAVu-Vudy< I -rj^AVu-Vu + 2{u- kfAV-q -V-qdy. 

Subtracting the first term on the right and using the elhpticity of A gives 

^[ Ti^\Vu\^dy<\ [ {u-kf\Vr^\''dy. 

Recalhng the estimate IVt^I < CQ/{as) and that 77 = 1 on Bs-as{y), we conclude 
that 

/ \Vu\^dy<j^\Ak,s\^Mu-kf- 
Thus u is in the space B{^r{x), 7) with 7 = ACq/X^. □ 

Lemma 3.3. Let x E dfl and assume that < r < tq. If u E B{^r{x),'y) and 
k is as in the definition of this space, then there exists 9i — ^1(71,7) .such that if 
\Ak,r\ < Oir" and H = sup^^ ^(w — fc) > 0, then |^fe+ir/2,r/2| = and hence 

sup (u - fc) < H/2. 



Before giving the proof of Lemma [3.31 we need to give two versions of the Sobolev- 
Poincare inequality. 

Lemma 3.4. Let x E and let r satisfy < r < tq. If u E W'^'^{'i/r{x)) and 
i > k, then 



(3.4) (^-fc)|A,.,|^-V"<C / \\/u\dy. 

I*r-(a;) \ Ak,r\ JAk,r\At,^ 

If B.r/2{x) n D 7^ 0, It e W'^^^{^r{x)), and k > 0, then 

(3.5) (i - k)\Ae^r\^-^/" < C f \Vu\dy. 

Proof. In each case, the estimate follows by applying a Sobolev inequality to the 
function 

{u{y) ~ k)+, u{y) < e 

i-k, u{y)>e 

The estimate p.5p uses our assumption that D satisfies (|2.2p . See Section 3 of [26] , 
for example. □ 



Proof of LemmalKM Fix ^(^(2;) and let u E B{'^r{x),j). For /i = 0, 1, 2, . . . , set 
r,, = r/2 + r/(2''+i) and kh ^ k + H/2 - H/{2''+^), 
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where H is as in the statement of Lemma 13.31 We will use the notation ah = 

[rh - rh+i)/rh. Since u e 7), we have 

(3.6) / \Wu\^dy<j- sup (u^kh)^ <j^^H^\Ak„r J. 

We use the inequality p. 41) from Lemma [3^ to obtain that 

(3.7) {kh+i - fc/.)|A,,^„r.+, r-^/" <C [ \Vu\ dy, 

where we choose 9i small in order to obtain a uniform bound on the constant in 
the Sobolev inequality in (13. 4p . Now p.6p and \?>.7\ give that 

^l^^c.,„,.,+,r-^/" < C I \Vu\dy 




< C{ \^u\'dy] |A,„.,,,r/2 



Thus, we may conclude that 
(3.8) 



/ \A-kh+i,rh+i I \ ^ ^ ^1/202/1+4 I ^fcfe, '■hi 

According to Lemma 4.7 in the monograph of Ladyzhenskaya and Ural'tseva 
[l9l p. 66] if 9i is sufficiently small, then the recursion relation p.8p implies that 
lim;,^o, |Afc,,,J/r" = 0. □ 

We now give the main step in the proof of Holder continuity of solutions of 
the mixed problem. Before stating the result, we introduce the notation osc^m — 
sup^ u — vciIe u for the oscillation of a real- valued function it on a set E. 

Lemma 3.5. Let x e 917 and assume that < r < tq. Let u be a solution of Lu = 
in '^r{x) and suppose that u has zero data for the mixed problem on d'^rix) n dO.. 
Then there exists an integer s such that 

osc w < (1 — 2^^*^) osc u. 

Proof. Since oscu — osc(— w), it suffices to prove the Lemma for either u or —u. 
We will take advantage of this in the proof below. We define 

Mr — sup{u(a;) : x E ^'^(x)} 
rur — inf{M(a;) : x e ^'^(a;)} 

UJr — Mr — mr — OSC U 
Mr = \^{Mr^mr). 

In what follows, set uj — uj^r and 

Dt = ^J\/4^~w/2',2r \ ^M4r-t^/2'+i,2r, ^ — 1, 2, . . . , S, 

where s remains to be determined. 
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There are a few details that are different in the cases when 'i>2r{x) n D ^ and 
^2r{x) n D = and we will point out the differences when they arise. 
In the case when '^2r{x) n = 0, we may assume that 

(3.9) \AM,^,2r\ < ^l*2r(a;)|, 

for if the condition (13.91) fails, we may replace u by —u. We next use the inequality 
(133) with k = AUr - and i = A'hr - w/2*+i to conclude that 

(3.10) ^I^M4,.-c./2*+i,2.r-'/" <C \Vu\ dy. 
Now since Dt C ^A/4^-tj/2*.2r and u G 6(^'r.(x), 7), we have that 

^u\^ dy < ^\AM„.-u./2^ir\ sup (m-(M4.-w/2*))2<c(|-)' 
Thus, from p.lOp and the Cauchy-Schwarz inequality, we have that 

(3.11) (^) Va/..-^/2*+^2.|'-'/" < CIAI (1)^"-'- 
If we sum (|3.1ip from t — l,...,s — 3, we conclude that 

s-3 

(3.12) {s - 3)|^M..-./2=-,2.p-'/" < Cr""2 ^ lAl < Cor2"-2. 

t=i 

Choose s such that 

^ ^ ri/(2n-2) 

^ < 01, 



2 



s-3 



with 61 as in Lemma 

, ./os-2 „„j o- ^^^^^ 



Now let A: = M^r - a;/2^-2 and = sup^ m{u ~ k) = Mr - {IVhr - If 



H > 0, we may apply Lemma 13.31 to obtain 

sup It < A: + H/2 < Mir - a;/2''"^ + ]-{Mr - M^r - a;/2^"^). 

*r(a;) 2 

Simplifying, we find that 

Mr < Mir - 

This inequality also follows easily if iJ < 0. It is immediate to see that —niir > —rrir 
and if we recall that lu = 1^4^, we may conclude that 

Next we consider the case when ^2r(a;) n Z3 ^ 0. In this situation we use the 
freedom to replace u by —u to impose the condition that Mir > and, as a result, 
p.9p is not guaranteed to hold. 

Since Mir — 1^/2* > A/4,. > 0, we may use the Sobolev inequality p.Sp to obtain 

^I^M.„-./2S2.r-'/" < C |Vu|dy, 

where D[ — AM^^_^j^t 4^ \^M4^_^^2t+i This replaces p.lOp in the argument that 
leads to p. 121) . The rest of the argument goes throughout without change. □ 

Proof of Theorem \3.1\ The Theorem follows immediately from Lemma 13.51 □ 
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We now return to working with only the Laplacian as this will simplify the 
uniqueness argument below. We define a Green function with pole at x for the 
mixed problem to be a function G{x,-) G M^^^^(il) which satisfies a) G{x,-) G 
Wj^ {n \ Br{x)) for each r > and b) if is in C°°{il) and vanishes on D, then 

(3.13) I S/G{x,-) -S/cjidy ^ -(j){x). 

The Green function is unique. If there are two candidates for the Green function 
with pole at x, Gi{x, •) and G2{x, •), then u = Gi{x, ■) — G2{x, ■) will satisfy 

Vu-V(j)dy = 



n 

for all (f) G C°°(f2) which vanish on D. From Weyl's lemma, u is smooth in the 
interior of ft. Then the assumption that each Gi{x, •) lies in Wjj (fl \ Br{x)) for 
each r > implies that u is in W^'^(n) Then M is a weak solution of the mixed 
problem with zero data and hence u = 0. The properties of the Green function for 
the mixed problem that we will need in the sequel of this paper are summarized in 
the following Lemma. 

Lemma 3.6. Consider the mixed problem in a Lipschitz domain with D satisfying 
112. S^) . Then there exists a Green function for the mixed problem which satisfies: 1) 
IfG^iy) = G{x,y), then G^ e W^^{n \ Br{x)) for all r > 0, 2) AG^ = S^, the 
Dirac S-measure at x, 3) If fn G ^^^'^(cJSl), then the weak solution of the mixed 
problem with = can be represented by 

u{x) = ~{fN,Gx), 

4 ) The Green function is Holder continuous away from the pole and satisfies the 
estimates 

\G{x, y) - Gix, y')\ < ^ , \x - y\ > 2\y - y'\, 

C 

\G{x,y)\<- — ^, n>3, 

and 

\Gix,y)\<Gil+log{d/\x-y\)), n = 2. 
Above, the exponent (3 is as in Theorem lS.li 

We give a detailed proof of this Lemma when n — 2. The proof for n > 3 may 
be obtained by adapting this argument or by a straightforward adaptation of the 
arguments in Griiter and Widman |12) . 

Proof. We begin with a result of Kenig and Ni [T7] who established the existence 
of global fundamental solutions in two dimensions. Suppose that L = divAV is 
an elliptic operator in two dimensions with bounded, measurable, and symmet- 
ric coefficients. Then there exists a fundamental solution T(x, y) which satisfies 
T{x, •) e Wl-^iR' \ {x}), r{x, •) e Wl^'jR'), and we have 

/ ^Vr(a;, •) ■\/(l)dy:^ for aU € C~(R2). 
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Chanillo and Li [4, Corollary 1.1] observe that the free space fundamental solution 
lies in BMO{Il^) and thus if F = r(a;, y) dy, we have 

/ {r{x,y)~ffdy<C 
JbAz) 

where the constant depends only on the ellipticity constant for the operator L and 
the bounds for the coefficients. 

Next we recall a result of Dahlberg and Kenig [9l p. 447], if VL' = {[xi,X2) ■ X2 > 
il>{xi)} is the domain which lies above the graph of a single Lipschitz function, 
then there exists a Green function for the Neumann problem in 57' which can be 
constructed by the method by reflection. We briefly recall the construction of this 
Green function. Define R a reflection in dVl' by R{x) = {xi,24'{xi) — X2) and 
construct an operator L = divAV on so that L{u o R) = Q m. \ f2' if and 
only if Au = in fi' and L = A in Let F be the fundamental solution for this 
operator in R^ and define 

N{x,y)=T{x,y)+T{Rx.,y). 

We have that A/" is a fundamental solution in U! with zero Neumann data on dO! . 
More precisely, we have the weak formulation 

(3.14) / yN{x,-)-^(l)dy^-(j){x), for all (/) e C~(f)). 
Jvt' 

We give a detailed proof of the construction of the Green function for the mixed 
problem G(a;, •) when x is near the boundary and hence lies in some coordinate 
cylinder. When x is far from the boundary, the construction of the Green function 
is simpler and we omit the details. Fix x in our original domain J7 and assume 
that x lies in a coordinate cylinder Zr{xo) with Ziri^o) also a coordinate cylinder. 
Let s = min(dist(x, -D), r) and then construct a cutoff function rj which is one on 
Bg/2{x) a-iid zero outside Bs{x). We let Af be the Green function for the Neumann 
problem in a graph domain O' that satisfies f2 H Z^r = f^' H Z^r- Since F and hence 
Af lie in BMO{Ii^), we may choose a constant N so that 

(3.15) / {J^{x,y)-Affdy<C, 

where the bound C depends only on the Lipschitz constant M . 

We will look for the Green function for the mixed problem in the form 

G{x, y) = r]{y){N{x, y) ~ Af) + u{y). 

We will show that the function u lies in W^'^i fl) with a bound depending only on 
the Lipschitz constant Af and the constant in (|3.2p . 
If G is to be a Green function, we need u to satisfy 

/ {-qWAfix, •) + {Af{x, •) - N')Wr] + Vu) -Vcfydy^ -<P{x) 
Jn 

for all functions (f> which lie in C°°(f2) and vanish on D. We write rjWcj) - 
(pS/rj and get 

/ WAf{x, ■)-V{T](l))-(l)VAf{x, •)• V7/+(7V(x, ■)-Af)VT]-V(j> + Vu-Vcj)dy 
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Since A/" is a Green function for the Neumann problem, we may use p.l4[) to simplify 
the previous equation and obtain 

(3.16) / Vu-V0rfy= / 0VA/'(a;, •)-V7?-(A/'(a:, •)-A/')V?7-V(/)dy. 

Jn Jn 

Let -F(0) denote the right-hand side of (|3.16p . We claim that |F(0)| < C||0||^i,2^j^^ 
where the constant depends only on the constant M and the constant in the coerciv- 
ity estimate p.2p . From the claim and basic results about Hilbert spaces, it follows 
that there exists a solution u to (|3.16l) and this solution satisfies ||w||,^i.2^^j < C. 

We now turn to the proof of the claim. To estimate the first term of p.l6p . we 
begin with an application of the Cauchy-Schwarz inequality 



(j>WAf{x, •) • Wrjdy 



< C 




4>^dy 



Using the Caccioppoli inequality and that M is in BMO{T{?) (see (|3.15p ). we may 
conclude 

/ „ \ 1/2 



/ \VU{x,-)\'dy 

where the constant depends only on the Lipschitz constant for VL. Thanks to the 
choice of s, we may use that (j) vanishes on D and (|2.2p to obtain the Poincare 
inequality 



/ <f>^dy<C f 

J'S'Jx) 



'dy<C 

See Section 3 of Ott and Brown [26j for details. 
Thus we obtain the estimate 



iV^l^dy. 



Vry dy 



4>VN{x, 

The estimate for the other term, 

(7V(x, •)-A/')V77- V(/)dy 



< C||V0||L2(n), 

follows from the Cauchy Schwarz inequality since N{x^ ■) is in BMO(R?). 

Next we recall that if u is in 1/Fi'^(17), x e Jl, and r > 0, we may find a constant 
u such that we have the Poincare inequality 



/ 



[u-ufdy < C||Vu|'2 



L2(n)- 



In other words, u lies in BAIO{fl). Since we also have that Af is in BMO(R^), we 
may conclude that G lies in BMO{^). 

Now we turn to the estimates in part 3) and 4) of Lemma [5751 First, recall that 
if V is in BMO{^) and '^r{x) n '^2r{x') ^ 0, then we have 



-f V dy — f 



V dy 



< V 



'*2r(a:') 

where ||f ||* is the BMO norm of v. Using this, an iteration argument and the local 
boundedness estimate p.3p . we obtain the pointwise upper bound, 

\Gix,y)\<Cil + log{d/\x~y\)), 
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where d is the diameter of il and the constant C depends on fl. Next we show that 
the Green function is Holder continuous. Let w be a solution of = in 'i'r{x) 
with zero data for the mixed problem on dflHd'^rix). From the local boundedness 
result (13. 3p and the estimate for Holder continuity in Theorem 13. 1[ for any constant 
V we have 

Hy)~viy')\<C{\y~y'\/rrJ \v-v\dy, y,y' e ^r/2{x). 

Since G is in BMO{Vl), the Holder estimate for G in part 4) of Lemma follows 
by applying the above observation to G{x,-) on a ball centered at y with radius 
r=\x- y\. 

Next, we claim that if € W^^^'^'^ {dVl) and u is the weak solution of the 
mixed problem with Neumann data /at and zero Dirichlet data, then we have the 
representation formula in part 3) of Lemma 13.61 

(3.17) = -(/jv,G,). 

Here, Gx — G{x, ■) and (•, •) denotes the duality pairing between W]j^^^'^{dVL) and 
W^J'^''^{dVl). To begin the proof of (jXTTI . consider VG(x, •) • ^udy. Fix xe^, 
let r — dist(a:, dVL)/2, and let 77 be a cutoff function with 77 = 1 on Bj./2{x) and rj = 
outside Br{x). We may approximate riG{x, •) in W-^'^{fl) by a sequence of smooth 
functions and use that u is harmonic in Br{x) to conclude that J^W{riG{x, ■)) ■ 
Vudy — 0. Since {l — r])G{x, •) lies in W]^'^{Q), we may use that u is a weak solution 
of the mixed problem to conclude that /qV((1 — 'q)G{x,-)) ■ Vudy = {Jn^Gx}- 
Combining these observations gives 

(3.18) / VG{x,-)-Vudy^{fN,Gx). 

Jn 

We now reverse the roles of G and u. With 77 as above, write u = rju+ {1 — rj)u. 
As u is harmonic and hence smooth in the interior of fi, there exists a sequence 
of smooth, compactly supported functions which converge to r]U in W^'°°{^). We 
may use this sequence and (I3.13P to obtain VG • V^rju) dy = —u{x). As u lies in 
W]^^{n) we may find a sequence {uk} of smooth functions which vanish on D and 
which converge in W^'^{ft) to u. Approximating (1 — ri)u by (1 — r])uk and using 
(13.131) gives VG(a;, •) • V((l — ri)u) dy = 0. Combining these observations implies 

(3.19) / VG(x, ■)-\Iudy^ -u{x). 

Jn 

From ([XTS|) and we obtain ([XTT]) . 

□ 

The next two lemmas establish higher integrability of the gradient of weak solu- 
tions to the mixed problem. The proofs of these lemmas appear in Ott and Brown 
[261 Section 3] . The key ingredients of the proofs appearing in the aforementioned 
paper are Poincare inequalities, and these Poincare inequalities continue to hold 
true in the current setting due to the condition (|2.2p on _D. A similar estimate is 
obtained for the mixed problem by Groger PJ using the method of N. Meyers [12] . 
However, Groger's method requires more restrictions on the boundary between D 
and N. 
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Lemma 3.7. Let ft and D satisfy and 12. S^) . Let x ^ Q, and let r satisfy 

< r < ro . Let u be a weak solution of the mixed problem for a divergence form 
elliptic operator with zero Dirichlet data and Neumann data f^ G L'^{N). Then u 
satisfies the estimate 




Above, q = 2 if n = 2 and q =^ 2{n — l)/(n — 2) for n > 3. The constant C depends 
on M and n. 

Lemma 3.8. Let ft and D satisfy \1.S\) and Ii2.2\) . Let x e 57 and let r satisfy 
< r < tq. Let u be a weak solution of the mixed problem with zero Dirichlet 
data and Neumann data f^ G L'^(N) which is supported in N D Ar{x), with q as 
in Lemma \3. 7| Then there exists go = 9o {M, n) > 2 such that for t in the range 
2 < t < qo when n > 3 and t in the range 2 < t < qo when n — 2, u satisfies the 
estimate 




The constant above depends on M and n. 



4. Estimates for solutions with atomic data 

In this section we consider the mixed problem with Neumann data an atom for 
N and zero Dirichlet data. We estimate the decay of the solution of this mixed 
problem as we move away from the support of the atom by taking L^-norms of the 
solution in dyadic rings around the support of the atom. Thus, given a surface ball 
Ar(a;), we define — l!^2*'r{^) \ ^2>'-'^r{x) s-^d = ^'21=^(2;) \ ^2''-^r{.'^)- 

Theorem 4.1. Let VL and D satisfy U.2]) and Fix x e dft and let r satisfy 

< r < tq. Let u be a weak solution of the mixed problem U.l]) with Neumann 
data /at = a an atom for N which is supported in (x) and zero Dirichlet data. 
Let qQ > 2 be as in Lemma \ 3.8\ and let A satisfy U.3\) with e such that < e < 
(go — 2)/(go — !)• Then for 1 < p < qQ{{l — e) / (2 — e)) , the following estimates hold 

(4.1) ( f \Vu\Pda\ < Ca{Arix))-^^P\ 



and for k> A, 
(4.2) 



j \Vu\Pda^ < C2'f^''aiY, 



i/p 



Here, /3 is as in Lemma \3.(A and the constants in the estimates and |^.^[ ) 

depend on p and the global character of the domain. 
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In order to prove Theorem 14. II we need a series of lemmas, some of which require 
that we work in a subdomain of fl which is also contained in a coordinate cylinder. 
For X e dft and r satisfying < r < ro, let flr{x) = Z^ix) n fi, where Zr{x) 
is a coordinate cylinder as defined in Section [2] The sets fir (a;) are star-shaped 
Lipschitz domains and for this reason they are preferable to the sets 4'r(a;)- 

The following lemmas rely on a Whitney decomposition of dn\A. For simplicity, 
we use surface cubes rather than the surface balls used up to this point. A surface 
cube is the image of a cube in R"~^ under the mapping x' — >■ (a;', (f>(x')). Then we 
write dVl ~ AlJ(Uj(5j) where the collection of surface cubes {Qj} has the following 
three properties: 1) For each j, either Qj C I? or Qj C (N \ A), 2) There exist 
constants c' and c", with c" as small as we like, such that for each x £ Q and each 
j, c'5{x) < diam(gj) < c"d{x), 3) If r(Qj) ^ {x e n : dist(a;,arj) < diam(Qj)}, 
then provided that the constants in the previous condition are sufficiently small, 
the sets {T{Qj)} have bounded overlaps and thus 

Y,XnQ,)<C{n,M). 

3 

We now begin the series of lemmas. The first two lemmas give a local version 
of boundary regularity for the Dirichlet and Neumann problems for the Laplacian, 
and they require only that i7 be a Lipschitz domain. The proofs of the next two 
lemmas appear in the previous work of Ott and Brown [26t Section 4]. 

Lemma 4.2. Let VL he a Lipschitz domain, let x € d^, and assume that r satisfies 
< r < tq. Let u he a harmonic function in Vlnr{x). If^u £ L"^ {Q,4^r{x)) and 
du/dv e L'^{d9. n d^ir), then Vu G L'^{Ar{x)) and 

f {{yu):fda<c(j p da + - f \Vu\^dy]. 

Ja^{x) \Jdandn4.^{x) r Jo4,.(:e) / 

The constant C depends on M and the dimension n. 



Lemma 4.3. Let he a Lipschitz domain, let x e dQ and r such that < r < tq. 
Let u he a harmonic function in fl^rix). IfVu G L^(ri4r(a;)) and Vtu G L^{dfl f] 
dfl4r{x)), then Vu G L'^{Ar{x)) and 

[ iiVu)*fda<cl[ \Vtu\''da+- [ \Vu\^ dy] . 

JA^ix) \Jdnndn4r{x) t JQ.i,.{x) j 

The constant C depends on M and the dimension n. 

The next lemma also appears in Ott and Brown [551 Section 4]. For the sake of 
completeness we sketch the proof again in this paper. 

Lemma 4.4. Let Q, and D satisfy U.S^) . and Let u be a weak solution 

of the mixed problem with Neumann data fj^ G L^{N) and zero Dirichlet data. Let 
p G R, a; G d^, and < r < rg, and assume that for some ^ > 0, 5(x) < Ar. 
Then it follows that 

[ ((Vu):,)2 S'-Pda <C{ [ \fN\'S'-P da+ [ \Vu\^ S'^ dy] , 

for constants c and C which depend only on M , n, A and p. 
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Proof. When Ar{x) is close to A, using the Whitney decomposition constructed 
above and the estimates of Lemma 14.21 and 14.31 we have 



(4.3) / {{Vu):fda<C 





du 




1 


dv 


da + — [ \\7ufdy\ 


)2QjnN 





To finish the proof of the Lemma, we multiply ()4.3|) by rj'' , recall that rj d{x) 
for all X e T{Qj), sum on the Qj that intersect Ar{x), and use that the family 
{T{Qj)} has bounded overlaps. □ 

The next result is another reverse Holder inequality, this time at the boundary. 
While at first glance the result below may not resemble a reverse Holder inequality, 
in future applications in this paper /^r = or a constant. 



Theorem 4.5. Let and D satisfy and \l-4^ - Let qo > 2 be as in Lemma 

\3.8\ and let A satisfy IIL3\) with < e < (50 ~ 2)/(qo — 1)- Let u he the weak solution 
of the mixed problem with Neumann data fi^ G L'^{N) and zero Dirichlet data. Fix 
p such that 1 < p < qo{l — e)/(2 — e). For x G dfl and r satisfying < r < Vq, 



1.4) (-f 



1/p 



(4.4) \4 \Vu\Pda] <C\4 |Vu|dy + r(i-")/2|l/. 



rA'||L2(A2,(x)nAf) 



Above, the constant C depends on M, the dimension n, and p. 

Proof. Fix X G dO, and < r < rg. We claim that 
(4.5) 

(/ \VuY'd<j\ <c{-f |V7.|dy + r(i-")/2|l/^|U.(^3^^(,)n^) 



We win separate the proof into two cases: a) b[x) < 8r V 1 + A'P , and b) 5{x) > 
8r\/l + M^. Starting with case a), choose p satisfying 2— |(1 — e) — e < p < 
2 — ^(1 — e) — e (the assumption that < e < {qo — 2)/{qo — 1) ensures that this 
is a non-empty interval). Apply Holder's inequality with the exponents 2/p and 
2/(2 -p) to get 

\ i/p 
\Wu\Pdcr] 

A4,.(x) / 

1 1 — 1 

< (/ \\/u\'^S^-P da] If 6^P-^^^da] 
(4.6) < Cr("-i)(p-^)+'^ f / \\'u\'^S^-Pda] , 



and Lemma [2^ to ensure that the integral of (5(''^^)(J'/(^~p)) is finite. Next, we use 



where we have used that (p — l)p/{2 — p) < —(1 — e) or p > 2 — ^(1 — e) — e 
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Lemma 4.6 and our hypothesis that d{x) < SrvT+M^ to obtain 

\ 1/2 



< c 

(4.7) < C 



\Vu\'^5-Pdy 



*Sr(2;) 



1/2 



1/2 



1/2- 



\fN\H^-Pda 



\\/u\^S Pdy] +r"2'"||/jv|U2(jvnA8.(x)) 



To estimate the term (J^^ |Vup(5 PdyY^'^, choose q > 2 such that q < 
minjqo, 2n/(n — 1)} and apply Holder's inequality again with exponents q/2 and 
q/ {q — 2) to conclude 

\ 1/2 
iVul^S-Pdy] 



*8r(a:) 



< ( f iVul^dy^ (f S-P^dy 



1_1 
2 q 



We invoke Lemma 13.81 again, requiring that q e (2, q^), and get the following bound 

1/2 



(4.8// \Wu\H-Pdy] 

y^Srix) J 



< Cr— 



\Vu\dy- 



g(Tl-l) 

|/jv| " da 



Ai6r(3:)nAf 



By our choice of q, we can apply Holder's inequality with exponent 2n/{q{n — 1)) 
to the boundary term in (14.81) to obtain 



(4.9) 



1 



^n— 1 / 

' JAi6r(a;)nAf 



" da 



< C 



1/2 



^n — l I 

' J Al6r-ix)nN 



IfNl^da 



Combining equations (|i?71) . and we conclude that 



7 



^ ( / |Vu| dy + r(i-")/2||/^||^.(A,,,(.)niv) 



which gives the claim (|4.5I) . 

Now we prove the claim (|4.5p under the condition in case b). As with case a), 
we begin with an application of Holders inequality 

l/p / ^ 1/2 



'Air (a;) 



j \Jau^) 
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Use that S{x) > SrvT+Af^ and Lemma [2.11 to conclude that either Ag,r{x) C N 
or Asr{x) C D. Then we may appeal to Lemma or Lemma to obtain that 

|2 ^ [ I I -P_j2 I ^ I IV7„,|2 



[sJuYd<T<C\ \fNrda+- \Vu\Uy]. 

Air-ix) \JAsr{x)nN ^ J'i'Srix) J 

From this point we must distinguish between n — 2 and n > 3. First, let n > 3. 
Then Lemma [3.81 and Holder's inequality give that 



/ \Vu\^dy] <C 

'■i'Srix) I 




\Vu\dy+ I / UNf^da 



'*i6r(a;) \J Aier{x)r]N ) 

Combining the last three displayed inequalities, we obtain 
[\ \VuYda\ <Cr^\l |Vu|dy + r(i-")/2|l/^r|U2(A,,,(,)n^) 

\JAir(x) / \J^mr{x) . 



which leads immediately to the claim (14.51) . In the case n = 2, we need an ad- 
ditional application of Holder's inequality. Choose t satisfying 2 < t < q^. Then 
(i*i6 [x] dyY^^ < C{f^^^ 1^^^ |Vm|* dyY/^ and from this point we can now ap- 

ply Lemma 13.81 to get the average of the square of the Neumann data. This gives 
(SSI). 

Once (14. 5|) is established, an elementary covering argument leads to the desired 
estimate. □ 



The final result we require before proving Theorem 14.11 is an energy estimate. 
The proof appears in Ott and Brown 26, Section 4]. 

Lemma 4.6. Let u he a weak solution of the mixed problem with Neumann data 
/at. For n > 3, let G LP{N) with p = {2n — 2)/n. Then the following estimate 
holds 



\^u\'dy<C\\fN\\UNy 
If n = 2, let fpf G H^(N) and then the following estimate holds 

\Wu\^dy<C\\fN\\j,i^j,y 

In both cases, the constant C depends on the global character of^. 
We are now equipped to prove Theorem 14. II 

Proof of Theorem \4-.l\ Fix x G d^l and let r satisfy < r < rp. The first step is 
to obtain an estimate for the gradient of the solution u near the support of the 
atom. Estimate (|4.ip follows immediately from Theorem 14. 5[ Lemma and the 
normalization of the atom. 

The next step is to estimate /j,^ \VuY' da for fc > 4. We begin by proving that 
the solution u satisfies the upper bound 

Crl^ 

(4.10) <____, \x~y\>2r, 
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where /3 is as in Lemma [3761 To establish ()4.10p . we use the representation formula 
in part 3) of Lemma 13.61 and claim that there exists x in /S.r{x) such that 



If Ar(x) C N , then let a; = a; and use that a has mean value zero to obtain the 
estimate (|4.10l) . If Ar(x) n £> 7^ 0, then choose x & Dd Aj.(a;) and use that G(y, •) 
vanishes on D. Now estimate (j4.10p follows easily from the normalization of the 
atom and the estimates for the Green function in part 4) of Lemma 13.61 

The remainder of the proof of estimate (|4.2|) follows from Theorem 14.51 and 
estimate ()4.10p . The constant in the estimate will depend on p, M, the dimension 
n, and the collection of coordinate cylinders. □ 

Next we prove that the non-tangential maximal function of a weak solution lies 
in (9ri) when the Neumann data is given by an atom. We introduce the following 
notation to be used in the proof. Let 

Ct = {y■.t<5{y)<2t}C^^ and Ct ^ {y : t < 5{y) < 2t} C^ d^. 

Theorem 4.7. Let and D satisfy \ 1.2(1 and Let A satisfy lil.3\) with 

< e < (go - 2)/(go - 1) with qo > 2 from Lemma lTR If fN G H'^{N), then there 
exists a solution u of the L^ -mixed problem hl.l]) with Neumann data fpf and zero 
Dirichlet data. This solution satisfies 



The constant C above depends on the global character of the domain. 

Proof. To begin, let /jv = a be an atom for N and let u be the weak solution of 
the mixed problem with Neumann data /at and zero Dirichlet data. The H^{N) 
estimate will follow immediately from the estimate for an atom. 

We wish to establish a representation formula for the gradient of u in terms of 
the boundary values of u. Let x d Q and let j be an index ranging from 1 to n. 
The claim is that 



where S is the fundamental solution of the Laplacian. In the case that u is smooth 
up to the boundary, the formula follows from the divergence theorem. However, it 
will take more work to prove (j4.1ip when u is only a weak solution. 

Let 77 be a smooth function that is zero in a neighborhood of A and supported 
in a coordinate cylinder. Using the coordinates of the coordinate cylinder, let 
Uriy) — u{y + re„), where e„ is the unit vector in the n-th direction. Applying the 




\\{^ur\\man)<C\\fN\\m[N). 



(4.11) 
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divergence formula gives 




dE 

(4.12) + — (x - OVur • Vrydy. 

Since the cutoff function t] vanislies near A, we may use tfie truncated maximal 
function estimate in Lemma l4.4l to let r approach from above and we can conclude 
that the identity (|4.12p continues to hold with replaced by u. 

Our next step is to remove this restriction that t] must vanish on A. Towards this 
end, suppose that r/ is of the form ij — rjcjjt, where 0t = on the set Cf , 0t = 1 on 
VL \ Czt and |V0t(x)| < C /t. By our conditions on the dimension of A and Lemma 
12.51 we have that 



(4.13) \Ct\ < Cf 



According to Lemma [3.81 Vu lies in L''(il) for q < qo- Using Holder's inequality 
with q < go, our estimate for V0t, and ()4.13p . we obtain 

(4.14) < Ct(2-^)(i/9')-i ^ |vu|9dy^ 



where 1/q + 1/q' = 1. The last term will go to zero as t goes to zero, if Vu E 
LiiQ) and (2 - e){l/q') - 1 > or e < (g - 2)/(g - 1). Our assumption that 
e < (go ~ 2)/(go ~ 1) implies that we may find a q for which the right-hand side of 
(|4.14p vanishes as t tends to zero. 

The remaining terms in (|4.12p can be estimated in a similar manner, yielding 

lim ' "'-^ " * " * 



[ Oil ( 

- / V(0t77) • \/yE{x -■)- \/yE{x - •) • Vii- 

Jn ciyj 



t^o+ Jn dy, dyj 

du „ s „ dr] 



-^{x - ■)Wu ■ W{(j)tii) dy 



.-y^v,.v.H(.-Oa;--v.E(.-.,.v„^ 

+ -;; — (a; — •)Vu • Vridy. 

Thus we obtain (j4.12l) with Ut replaced by u and rj is not required to vanish on A. 
Choose a partition of unity which consists of functions that are either supported in 
a coordinate cylinder, or whose support does not intersect the boundary of VL. As 
•q runs over this partition, the sum gives the representation formula for Vm given in 
(j4.1ip . From Theorem 14. II we have Vu G LP{dfl), and by the theorem of Coifman, 
Mcintosh and Meyer ^ it follows that (Vm)* e LP{dn) and thus (Vw)* G L^{dn) 
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since 51 is a bounded domain. A few more steps will provide us with the desired 
estimate, ||(VM)*|jLi(an) < C. 

Since (Vw)* lies in LP{dn), we can apply the divergence theorem twice more and 
obtain the following identities 



Using these identities and the estimates for Vw established in Theorem 14.11 we 
can conclude that the two integrands above are molecules on the boundary of the 
domain, and hence it follows from the representation formula ()4.1ip that (Vm)* lies 
in L^{dft) and satisfies the estimate 



The estimate for solutions with Neumann data in H^{N) follows easily from the 



We now turn to establishing uniqueness of solutions to the mixed problem. We 
rely on the results of the previous section and uniqueness of the regularity problem 
due to Dahlberg and Kenig [9] (also, see the work of D. Mitrea [23, Corollary 4.2] 
for the result in two dimensions). More specifically, we prove that if u solves (ll.ip 
with zero data and (Vu)* G L^{dft), then u also solves the regularity problem with 
zero data and hence u = 0. 

Theorem 5.1. If u is a solution of the L^-mixed problem (1.1) with fo — o.'nd 
fN = 0, then u = 0. 

The proof of this theorem closely follows the proof of uniqueness in the paper of 
Ott and Brown |26| Section 5]. We outline the main steps of the argument here 
and will omit many of the technical details. 

Recall the following approximation scheme of G. Verchota (j33l Appendix A] 
and [34j Theorem 1.12]). Given a Lipschitz domain fi, there exists a family of 
smooth domains {fi/t} with Cl^ C 51, k — 1,2,..., and a family of bi-Lipschitz 
homeomorphisms A^ : dVL^ — > dVl. 

The following lemma is proved via the method of Verchota. The argument 
is sketched in Ott and Brown [26l Section 5]. The proof uses generalized Riesz 
transforms and also relies on the Holder continuity of the Green function established 
in Lemma 13.61 We use a similar argument in Section [71 

Lemma 5.2. Let {51^} he a family of smooth domains. Let u G W^'^{dnk) for 
k ^ 1,2, . . .. If w is a weak solution of the mixed problem in flk with Neumann 
data an atom for N and zero Dirichlet data, then 




\\{^urh^9n)<C. 



result for solutions with atomic data. 



□ 



5. Uniqueness 




Next we introduce a Poincare inequality that will be employed below. Let = 
{x £ N : 6{x) < e}. We show that there is a constant C such that for u G W^'^{dil), 
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with u = a.e. on D, 

(5.1) - [ \u\da<C [ \Vu\da. 

To prove the inequahty above, let x ^ and let Qx.t denote the surface cube 
centered at x with side length e. If x S A^e, then Qx^2<i n £> ^ 0, thus we may use 
(|2.2p to obtain the Poincare inequality 

\u\d(j <Ce f \\7u\da. 

By the Besicovitch covering lemma, we can find a finite cover of by cubes 
Qxi,4e^ Qx2A^^ ■ ■ • 5 QxmA^ such that the cubes have bounded overlaps, 

m 
i=l 

Then we have 



/ \u\d(j < V / \u\da 

JN, j = l"'Qxi.4e 

< C^e / \Wu\da 



< Ce I \Wu\d<j. 

JNca 

Thus wc obtain (|5.ip which we will use to prove the following approximation 
lemma. Let T denote the collection of functions defined in that are Lipschitz 
and compactly supported in dU, \ D. The next lemma shows that a function u S 
W^''^{dVl) which satisfies u — Q a.e. on D can be approximated in the W^''^{d^)- 
norm by functions in T. This density is fairly easy for the domains considered in 
Ott and Brown |26| . but requires more work under the present assumptions. 

Lemma 5.3. Let Q. and D satisfy IJ.g)) and i2.^) . Suppose that u G W^'^{d^l) and 
u — a.e. on D. Then u can be approximated in W^'^{dQ.) by functions from T. 

Proof. Let u e W^'^{d^l) and suppose that u — a.e. on D. Fix e > and 
let 7]^ be a smooth function which is 1 if f > 2e and if t < e, and satisfies 
|V77(:(a;)| < C/e. If m S W^'^{dil) and vanishes a.e. on D, then we have that 
ite(x) — rif{6{x))u{x) is zero in a neighborhood of D U A. Then it follows that 
VtWe — VfU = uWtrie{S{-)) + (?7e((5(-)) — l)Vfit. From the dominated convergence 
theorem, 

lim / \\7tu\\r],{S{-))-l\da^0. 

Since Vr/e(5(-)) < C/e, we may use the Poincare inequality (|5.ip . which requires 
|), and the dominated convergence theorem to conclude that 



lim / \u\/tvAS{-))\da = 0. 



Thus, we have that u G W^'^{dQ,) may be approximated by a function that is 
supported in dfl \ D. By a standard regularization argument, we can approximate 
We by functions that are in T. □ 
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Proof of Theorem 15. il Let u be a solution of the mixed problem (jl.ip with /at = 
and /d = 0. We wish to show that u — 0. Fix an atom a for N and let w be a 
solution of the mixed problem with f^ — a and /d = as constructed in Theorem 
4.11 Our goal is to show that 



(5.2) / uada = 



N 

In turn, this will imply that u is zero on dfl, and by appealing to the uniqueness of 
the regularity problem proved by Dahlberg and Kenig [9 or D. Mitrea [53] in two 
dimensions, we can conclude that u — in n. 

To prove (|5.2p . we apply Green's second identity in one of the smooth approxi- 
mating domains from Verchota's construction and obtain 

(5.3) / da ~ [ u^^da, fc = l,2, 

We have that (Vu)* e L^{dVt) and Lemma [3751 implies that w is Holder continuous 
and hence bounded. Further, w ~ Q on D and ^ = on A^. Hence by the 
dominated convergence theorem, 

(5.4) lim / w^da = Q. 
Thus, we can prove our claim by showing that 

(5.5) lim I u^— da ~ f uada. 

Note that the existence of the limit in (15.51) follows from (|5.3p and (15.41) . Now by 
repeating the argument used to prove Lemma 5.7 in the work of Ott and Brown 
[26] . we can find a sequence {Uj} of Lipschitz functions defined in f2 such that 
Ujlo = and 

(5.6) lim \\u- UjWw^.i^gn,) < l/j- 

The argument outlined above uses the density result in Lemma |5.3[ Now we have 
I / uada— lim / u——da\ < | / uada— lim / Uj——da\ 



(5.7) 



lim sup I / {u — Uj)-:r— da\. 

k^oo JdVlk OV 



Since (Vw)* G L^{dVl) and Uj is bounded, we may take the limit of the first term 
on the right of ([ST]) . This yields 



/ uada— lim / u-—da\<\ f {u — Uj)ada\ + C/ j < C/ j. 
Jan '^^°°JdQ.k ov Jq^ 



The second term on the right of ()5.7|) is bounded by C^/j by Lemma 15.21 and 
(|5.6|) . Since j is arbitrary, we have obtained (|5.5|) and the proof of the Theorem is 
complete. □ 
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6. LP RESULT 

In this section, we use the existence of solutions of the mixed problem with 
data from Hardy spaces established in Section 2] to prove L^-estimates for the 
mixed problem. Our strategy is to first recall the reverse Holder inequality at the 
boundary which was proved in Theorem 14.51 With this estimate in hand, we then 
apply the method developed by Shen [28] and adapted by Ott and Brown [26] to 
obtain the L^'-estimate. 

The following lemma is a local estimate that is a consequence of Theorem 14.51 
In this Lemma we use the truncated non-tangential maximal function defined in 
Section d 

Lemma 6.1. Let and D satisfy assumptions U.S^) and ( [_?.^[ ). Let > 2 be as 
in Lemma \3.8\ and let A satisfy liL3\) with < e < {qo ~ 2)/ (go ~ !)• Let u be the 
weak solution of the mixed problem with f^ G L'^{N) and zero Dirichlet data. Let 
a; g O and < r < tq. Then for 1 < p < qo{{l ~ e)/(2 — e)) the following local 
estimate holds 

lA^ix) J \J-f2rix) J 

The constant c — 1/16 and C depends on M and n. 

Proof. Let x € ^ and r satisfy < r < tq. Theorem 14.51 provides an estimate for 
the LP-norm of Vm in a surface ball Ar{x). To obtain the estimate for the non- 
tangential maximal function, choose a cut-off function rj which is one on B^rix) 
and supported in B^rix). Let z G Br{x). By repeating the argument used to prove 
()4.1ip in the proof of Theorem I4.7| we can show that the gradient of the weak 
solution u may be represented as 

dzj Jgfj dv dyj dy, dv 

as 

-I-V77 • Vm— — iz — •) dy. 

From this representation and the theorem of Coifman, Mcintosh and Meyer [5] on 
the boundedness of the Cauchy integral, we get 



(6.1) {-f {\/uy/d<^ <C 
\Jaax) I 



\Wu\dy 




Prom estimate (I6.ip . Theorem 14.51 and a covering argument, we obtain the Theo- 
rem. □ 

Next, we outline the argument developed by Shen [5H] that we employ to obtain 
LP-estimates in this section and weighted i^-estimates in the next section. Shen's 
argument is adapted from work of Peral and Caffarelli [3] . It depends on a Calderon- 
Zygmund decomposition of the boundary and thus we will use surface cubes in this 
section rather than surface balls Ar{x). Before giving Shen's result, recall that a 
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locally integrable function w is an Ap(da) weight, 1 < p < oo, provided that 



(6.2) 



w da 



w-P'/P da 



p/p' 



< A < oo, 



a(A) \aiA) 

for all surface balls A C dfl centered on dil. Define Aao{da) — UpAp{da). 

Let Qo be a surface cube and let F be defined on AQq. Let the exponents p,q 
satisfy 1 < p < q. Assume that for each Q C Qo, we may find functions Fq and 
Rq, defined in 2Q, satisfying 

(6.3) \F\<\Fq\ + \Rq\, 



(6.4) 



(6.5) 



/ 

J2Q 



Fqlda < C 



J4C 



\f\'da 



i/p 



/ \Rq\ 

J2Q 



1/9 



\F\da- 



4Q 



da 



i/p- 



Going further, assume that /i is a weight in At{da) and that 



(6.6) 



i<t< 



where 9 depends on M. Under the assumptions (|6.3p - (|6.6p . for s in the interval 
{p, 0q), we have 



(6^7) 



/ 



F\' fida 



l/s 



< C 



\F\da- 



|/|>da 



4Qo 



l/s 



The constant in the estimate above depends on the Lipschitz constant of 17 and the 
constants in the estimates (|6.4p - (|6.6p . The argument to obtain (|6.7p is essentially 
the same as in Shen [551 Theorems 3.2, 3.4]. Ott and Brown [521 Section 7] rework 
Shen's argument to apply to the current situation, where our starting point is a 
result in a Hardy space rather than in an i^-space. 

Let AQo be a surface cube with sidelength comparable to tq. Let u be a solution 
of the mixed problem with Neumann data / £ LP[N) and zero Dirichlet data. 
Since LP{N) C H^{N), we know by Theorem 14.11 that a solution u exists with 
(Vm)* e L^{di}). Let F — (Vu)*. Now given a cube Q C Qo with diameter r, 
we define Fq and Rq as follows. Let /4Q = if 4(3 n I? 7^ and /4Q = / da 

if 4Q C N. Set g — XiQif ~ fio) and h = f — g. By construction, g and h are 
both elements of H^(N) and thus we may solve the mixed problem with Dirichlet 
data zero and Neumann data g or h. Let v solve the mixed problem with Neumann 
data equal to g and let w solve the mixed problem with Neumann data h. By 
our uniqueness result Theorem 15.11 we have that u — v + w. Let Fq — (Vf)* and 
Rq = (Vw)*. It follows immediately that (|6.3I) holds. The procedures for obtaining 
estimates (|6.4p and (|6.5p are straightforward and were worked out in detail in Ott 
and Brown [26 . 

With l|6.3p - (|6.5p established, we obtain (|6.7p for s in the interval {p,6q) and 
q < go((l — £)/(2 — e)). From here, we can now easily complete the LP-estimate. 

Proof of Theorem \l.l\ To prove part b), we use Dahlberg and Kenig's result [HI 
Theorem 4.3] for the regularity problem in Hardy spaces, or D. Mitrea's result [23] 



THE MIXED PROBLEM WITH GENERAL DECOMPOSITIONS OF THE BOUNDARY 29 

in two dimensions, to reduce to the case where the Dirichlet data is zero. Then we 
take /at G H^{N) and use Theorem 14.71 to complete the proof. 

To prove part a), consider the mixed problem with zero Dirichlet data and Neu- 
mann data in U'{N). Since LP{N) is contained in the Hardy space H^{N), a 
solution u exists by part b). From the argument of CafFarelli and Peral, as adapted 
by Shen, and Ott and Brown, we have that u satisfies the estimate 

||(VM)*||LP(an) < C\\fN\\LP(N)- 

Uniqueness of solutions of the mixed problem follows from Theorem 15.11 □ 

7. Weighted result 

In this section we establish results for the mixed problem with data from weighted 
Sobolev spaces. Throughout this section we assume that and D satisfy conditions 
(lL2)l . (fOl) . and (lO)) . 

To begin, we consider the regularity problem when the data comes from a 
weighted Sobolev space. We will use the solution of the regularity problem to 
reduce the study of the mixed problem to the case when the Dirichlet data is zero. 
Our study of the regularity problem contained here is a small extension of work 
of Shen [57] who studied the regularity problem with data in weighted L^-Sobolev 
spaces. Shen's work is in turn an extension of a method used by Verchota [33] to 
study the (unweighted) regularity problem in Lipschitz domains. This method is 
also developed in a recent article of Kilty and Shen 18 that studies the relationship 
between the regularity problem and the Dirichlet problem for elliptic systems. We 
choose to repeat well-known arguments for several reasons. Kilty and Shen do not 
give weighted estimates and there is a small mistake in [27]. The weight defined in 
equation (7.29) on page 2868 of [27 may not be a doubling weight and hence may 
not be in any Ap class. 

The heart of the matter is Lemma 17.31 below . which estimates the normal deriv- 
ative of a harmonic function in terms of its boundary values. Building toward this 
result, we begin by recalling that Verchota's result for the regularity problem with 
data in unweighted Sobolev spaces depends on a duality argument and the solu- 
tion of the Dirichlet problem with data in weighted i^-spaces. Thus, our starting 
point will be the following result of Dahlberg [7] regarding the weighted Dirichlet 
problem. 

In the results that follow, constants have the dependencies given in Section 2. 
In addition, the constant may depend on the weight through the Ap-constant and 
the exponents appearing in the assumptions on the weights. 

We begin by recalling some well-known results about the LP{fi (icr)-Dirichlet prob- 
lem. In this problem, given / on the boundary we look for a function u which 
satisfies 

{Am = 0, in rj 

u = f, on dn 

u* e LP{pda). 

Theorem 7.1. There exists an exponent sq < 2 such that if fJ. € Ar{da), r > 1, 
p > rsQ, and f G L^lpda), then the (p, da) -Dirichlet has a unique solution. 

Proof. Dahlberg [7] has shown that there exists an exponent to > 2 such that the 
harmonic measure lies in the reverse Holder class Bt{da) for t < to, meaning that 
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for each t < to, there is a constant Ct such that 
(7.2) uj^da^ - / 

for any surface ball A centered on dfl. Here, ut denotes the density with respect to 
surface measure of harmonic measure at some convenient point in fl. The exponent 
So will be the dual exponent to to, i.e., I/sq + 1/to = 1. 

Let / £ L^{da) n LP{^da), and let u be the solution of the L^((icr)-Dirichlet 
problem. From Hunt and Wheeden [13] (see also Jerison and Kenig [14|), we know 
that 

u*{x) <CM^f{x), 

where M^^f is the Hardy-Littlewood maximal function with respect to harmonic 
measure given by 



M^f{x) = sup— [ \f\ujda. 

r>0 Uj{Ar(x)) JaAx) 



r>o w(Ar(a;)) Ja,{x) 

Since w G Bt{da) for t <to,^e have 

M^f{x) < CsM{\fYY''{x), so < s < oo, 1/so + lAo - 1. 

Next, set s = p/r and note that our assumption p > rso implies that s > so- 
Since fi G Ar{da), it follows that with this choice of s we have the estimate 

\\M{\fn'/^\\LH^d.)<C\\fhn^a.), 

which implies that 

(7.3) IklUnM.) <C|1/I1lp(m.)- 

By a standard limiting argument, we may therefore construct solutions u £ LP{p da) 



to the Dirichlet problem (17. ip which satisfy the estimate (|7.3p . 

Finally, to establish uniqueness of solutions of (17. ip . observe that LP{fida) C 
L^'ida) when s — p/r. Since s > so, we may use the uniqueness result for the 
L''(d(T)-Dirichlet problem to conclude that if u is harmonic in J7, u* € LP{fida), 
and u has non-tangential limits of a.e. on dfl, then u — 0. □ 

The next theorem establishes solvability of the regularity problem when the 
boundary data lies in a weighted Sobolev space. Given a function / on the bound- 
ary, the LP {^J.~P Ip (i(T)-regularity problem is the problem of finding a function u 
which satisfies 

Au = 0, in 
u — f, on dfl 

(Vw)* e LP{fi-P'/P da). 



Theorem 7.2. Let so be as in Theorem \7.1\ and let fi G Ar{da), r > 1. If oo > 
p > rso and / lies in W^'P {fJ-^P ^Pda), then there exists a unique solution of the 
LP ip^P ^P da) -regularity problem with data f which satisfies 

{{Vu)*)P' p-P'IP da <C [ {\\7tu\P' + \u\P')p,-P'/Pda. 
an Jan 

In the following statement, let Vt denote the tangential gradient at the boundary 
(see Section 2 for the definition). 
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Lemma 7.3. Let /i G Ar{da) with r > 1 and suppose that oo > p > rso, where sq 
is as in Theorem \7. 1\ If u is a harmonic function with (Vu)* G L^{da), then 
du / \ 

\\'g^\\Lp'{fj.-p'/Pda) - ^ (^ll^t"llLp'(M-p'/PdCT) + 1 1 "I I Lp' (p-p' /p dcr) J ' 

To prove this lemma, we begin by defining local Riesz transforms. Fix a coordi- 
nate cylinder Zr{x), r < ro, for dfl such that Z2r{x) is also a coordinate cylinder. 
Hereafter in this section we will use Zr to denote Zr{x). Let 77 be a smooth cutoff 
function such that 77 = 1 on and rj — outside ^3^/2- Let be a harmonic 
function. Using coordinates (a;', a;„) G R"^^ x R, for i = 1, . . . , n, the local Riesz 
transforms are given by 

f°° d 

Viix) ^ - —{rjv){x',t)dt, 

Jx„ OXi 

for X G Z2r n fi, and Vi{x) = in i7 \ Z2n i — 1, . . . , n. Straightforward calculations 
give that 

dvi dvj 



dxj dxi ' 



(7.5) ^^(x) - / Vr^{x',t)-Vv{x',t)+v{x',t)^T^{x',t)dt 

Avi{x) = -/ —{\7r]{x',t)-\/v{x',t) + v{x',t)AT]{x',t))dt. 

Jx„ OXi 

Lemma 7.4. Fix x G Sfi and < r < rp. iei &e t/ie local Riesz 

transforms of a harmonic function v in a coordinate cylinder Z^2r o,nd suppose that 
fi lies in Aao{da). Then for p < 00, the following estimate holds 

I {vl,rfida<C{[ {virfida + fi{Asr{x))snp\v\P). 

Above, K C Z^2r is a compact subset of £7 and the cone opening for the non- 
tangential maximal function on the right is larger than the cone opening for the 
non-tangential maximal function on the left. 

Proof. The proof uses a truncated square function which we define by 
Sr{v){x) = (^J^ ^ ^\Vv{y)\'\x ~ y\'~- dy^ , 

where Tr{x) is a truncated cone as defined in Section [21 Let w be a given harmonic 
function and then let Vi be one of the local Riesz transforms of v defined in a 
coordinate cylinder Z^2r- Write Vi — v[ + where v[ is harmonic in i7 n ZiQr and 
v'/ — "E.* (AwiXZierno)- We observe that Avi is bounded in ZiQr and 

C 

sup \Avi\ < — sup \v\, 

where K is the compact set 

K = {{y', y„) : \y' - x'\ < 32r, (1 + M)r < y„ - x„ < (1 + M)32r}. 
Thus, we have 

(7.6) sup Kl +r|Vi;f I < Csup|w|. 

ZiQj. K 
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With these preliminaries, we can now give the main estimate 



C{[ {vl/r^ida + ^i{Ar{x))snp\v\n 




< 



C{ f S2r{v'i)'' f^da + ii{Arix))snp\v\P) 




< 



C{ [ S4rK)P fida + Ai(A,(a;)) sup \v\p) 




< 




The first inequahty follows from ()7.6p . the second is a local version of a theorem 
of Dahlberg 8j Theorem 1], the third inequality follows from a pointwise estimate 
which may be found in Stein |30l p. 213-214], and finally the fourth inequality 
follows from Dahlberg's result and (|7.6p . Note that in each of the inequalities 
above, the cone opening for the object of the left side must be smaller than the 
cone opening for the object on the right side. Our notation is already elaborate 
and thus, we choose to suppress this dependence. Once we recall that w„ = riv, 
then the estimate of the Lemma follows. □ 



We now are ready to present a proof of Lemma 17.31 



Proof of Lemma \7.3\ Let u be a solution of the (rfcr)-regularity problem with 
data / e W^''^{da). Since we may solve the ((i(T)-rcgularity problem [l^, we may 
assume that / is supported in a surface ball Ar{x), x £ dfl, and that Z32r = Z32r{x) 
is a coordinate cylinder. Suppose that dfl is given as the graph of (j> in ^32^. We 
would like to show that 



Toward this end, choose g G W'^''^{da) and let v be the solution of the Dirichlet 
problem with data g. We observe that 



Since v — g on dfl, u and v are harmonic, and (Vu)* + (Vw)* € Lp'{da), we may 
use Green's second identity to conclude that 



du 



< '^\W\\w^:v' /V da)- 



'Vp da) 
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As u = on dn\Zr, we may use that Vn — rjv = v on Zr, (|7.4p . (|7.5|) . some algebra, 
and integration by parts to obtain 

dv f dv„ 

u—— da — / u— — da 



dn Jan 



n-1 



i—i 



n— 1 ^ 

a 



u{x', (t){x')){F{x' ,(t){x')) + ^ —Vi{x', (t){x'))) dx 

2 — 1 * 



n-1 



where = (1 + | V(/)p)^^/^(ei + (pxi^n) is a tangential vector, Vi are the local Riesz 
transforms, and F — ff^ ^^e right-hand side of (|7.5p . Thus from Lemma 

41 and Theorem 17. 1[ it follows that 



f dv 

\ U—da\ < \\u\\Li{da)\\F\\L'^(da) + \\v*\\LPipda)\\'^tu\\Lp'(p.~p'/Pda) 

Jon 

— ^\\^\\w^'P'(tJ.-p'/p da) \ \9\\LP{iidcr)- 

We give the details for the estimate for the term ||w||ii((ic,)||F||ioo(^o.). By a Poincare 
inequality and Holder's inequality, we obtain 

\\u\\L^da) < Cr\\Wtu\\Li{da} < Cr || Vj u|| ^p' (^-p'/Pd^r) m( (a;) ) . 
Recall that F is the right-hand side of ()7.5p . Then it follows that 



WFh^^da) < Cr-^ sup\v\ < Cr-^ti{Ar{x))-^/P\\v*\\LP(j,d.). 

K 

With these inequalities, the stated estimate follows. □ 

Before proving Theorem 17.21 we require one more standard lemma. 
Lemma 7.5. // G Ar{da), 1 <r < oo, Au — 0, and (Vu)* G L^{da), then 

[ {{Vnyy nda<C [ iVu^fida. 

Proof. Let S be the fundamental solution for the Laplacian. We assume that 
(Vu)* G L?{da). Under these conditions, it is easy to establish the representa- 
tion formula 

du , . f dE du du dE sdu 

d^^'^ ^ " Jon - -^^^'w, ' % ^ % " % 

The Lemma now follows from standard estimates on singular integral operators 
on Lipschitz surfaces jj] and weighted estimates for Calderon-Zygmund operators 
0. □ 

We are now prepared to prove Theorem 17.21 
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Proof of Theorem\T^ Let / G T4^1'P' (^"p'/p da) D W^'^{da) and let u be the solu- 
tion of the ((i(T)-regularity problem with data /. From Lemma 17.31 and Lemma 
17.51 we conclude that 

iiVuff' ii-P'/P da<C [ i\Vtu\P' + \u\P')fi-P'/P da. 
en Jdn 

Now a limiting argument gives existence of solutions of the regularity problem with 

boundary data / e W^-p' {ij-p' /p da). Since Lp' {fi-P /p da) C L^{da), uniqueness 

follows from Dahlberg and Kenig [5]. □ 



Theorem 7.6. Let to be as in the reverse Holder inequality and sq — io/(^o~ 

1). Suppose that /i e Arida), where p > rsg. Set 

p-1 1 r-1 p-r 

a — and tl — — = 1 — -. 

r — 1 a p — 1 p ^ 1 

Assume that 1 < p' < 6{qa/2). Consider the mixed problem with Dirichlet data 
fo € W^'P {l^^P Ip da) and Neumann data /^v S LP {ijl~p Ip da). Then the following 
estimate for the solution u of the mixed problem holds 

(7.7) II (Vw)* II ^p'(^-p'/p do-) ^ ^ (IIInIIlp' (fi-p' /P da) + WfoWw^-p' (fi-p'/P da-)^ ■ 

Furthermore, there is only one solution satisfying (Vu)* G LP {pi~P ^p da). 

Proof. In our application, we are given /i G Ar{da) and p > rsg- We will apply 
(|6.7p with /i^P /P G Ap' [da) and we observe that with a — we have 

{fi-P'/P)°'da = / 

J A 

< [ I fida] a{Ay'' 

-l)/(r-l) 



A 



< 

for any surface ball A. Therefore 



U pi-P'IPda\ a(A)'-'a(A)-P'/('-i\ 



(7.8) J fi'-''P'°'^/Pda < (J fi-P'/Pda^ , 

where we have used the Ar condition for /i, 



and Holder's inequality, 

/ I- \v/(r-l) f f N (p-l)/(r-l) 

a{A)^/^'-^') < U fidaj Ufi-P'/Pdaj 

From the inequality (17. 8|) we have ^~p /p G Ba{da), and Holder's inequality 
implies 

^iiE)\ ^ / a{E)_y ^ Q^^_}__p-r 



^^(A) J - Vct(a) y ' a 1' 

Thus, we have the conditions needed to obtain the conclusion (16.71). 
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To summarize, if 1 < p' < 6qo/2, with go as in Lemma l3.8( and if /i G Ar{d(7), 
p > rso, with So as in Theorem lT-H then we have the solution of the mixed problem 
satisfies 

ll(V")*llLP(p-p'/Pdcr) ^ (\\fN\\Lp'(^^-p'/Pda) + WfoWw^.p'if^-p'/Pda)) ■ 

Since {fi"^ Z^) C L^{da), uniqueness for solutions of the mixed problem with 
(Vu)* S LP {fi'P ^P da) follows from the problem follows from part a) of Theorem 
1.1. 

□ 
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